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\q [ I sketch how long wavelength modes of the pion field can be ampli- 

fied during the QCD phase transition. If nature had been kinder, and 
had made the pion mass significantly less than the critical temperature 
00 ! for the transition, then this phenomenon would have characterized the 

■ transition in thermal equilibrium. Instead, these long wavelength os- 

. cillations of the orientation of the chiral condensate can only arise out 

| of equilibrium. There is a simple non-equilibrium mechanism, plausibly 

operational during heavy ion collisions, which naturally amplifies these 
oscillations. The characteristic signature of this phenomenon is large 



fluctuations in the ratio of the number of neutral pions to the total 
number of pions in regions of momentum space, that is in phase space 
in a detector. Detection in a heavy ion collision would imply an out of 
c~| , equilbrium chiral transition. 



Shortly after the discovery that QCD is asymptotically free |T| and that 
therefore quarks are weakly interacting at short distances, Collins and Perry B 
noted that this means that at temperatures T 3> Aq CD , the theory describes 
a world of weakly interacting quarks and gluons very different from the famil- 
iar hadronic world.0 There are at least two qualitative differences between a 
plasma at T ^> Aq CD and at T ~ 0. First, at low temperatures one has a 
plasma of hadrons, while at high temperatures the quarks and gluons are de- 
confined. Second, at low temperatures a qq condensate spontaneously breaks 

1 Talk given at the International Workshop on QCD Phase Transitions held in January, 
1997 in Hirschegg, Austria. 

2 In this talk I will limit myself to discussing the physics of a QCD plasma with zero 
baryon number density. This should be a good approximation in the central rapidity regions 
of heavy ion collisions at RHIC energies and higher, where the number of pions per event is 
expected to be about two orders of magnitude larger than the baryon number per event. 
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chiral symmetry while at high temperature the interactions among quarks and 
anti-quarks are weak, no such condensate exists, and chiral symmetry is man- 
ifest. 

In 577(3) gauge theory with no quarks, there are analytic arguments]^] con- 
firmed by lattice simulation Q that the deconfinement transition is first order. 
While this conclusion presumably remains valid if all quarks are much heavier 
than ~ Aqcd; no deconfinement order parameter is known for SU(N) gauge 
theory including dynamical quarks, and simulations with two light quarks || |6| 
find that deconfinement occurs via a smooth crossover. 

The presence of light quarks raises the possibility of a phase transition 
associated with the vanishing of the chiral order parameter. The Lagrangian 
for QCD with two massless quarks has a global SU{2) L x SU(2) R symmetry 
which, at low temperatures, is spontaneously broken to SU (2) L+R by a nonzero 
expectation value for the chiral order parameter, which can be written in terms 
of four real scalar fields (cr, 7?) and the Pauli matrices according to 



In fact, the order parameter can be written as a four component scalar field 



and the SU{2) L x SU{2) R transformations are simply 0(4) rotations. At low 
temperatures, (<ft) is nonzero. This picks a direction in 0(4) space defined to 
be the a direction, and spontaneously breaks the symmetry. The direction in 
which the expectation value points is defined as the a direction. Oscillations 
of in the 7? directions, that is oscillations of the orientation of the conden- 
sate, are massless Goldstone modes. As T is increased, fluctuates more and 
more wildly until above some critical temperature T c the fluctuations are large 
enough that (0) = and the 0(4) symmetry is restored. 

Since the order parameter appropriate for the chiral phase transition in 
QCD with two flavors of massless quarks has the symmetries of an 0(4) Heisen- 
berg magnet, and since this model has a second order phase transition,^ it is 
possible that the chiral phase transition is second order and is in the same 
universality class. |7|, [|, || At T = T c for a second order phase transition, the 
theory is at an infrared fixed point of the renormalization group and physics is 
scale invariant. This means that the order parameter fluctuates on all length 

3 This has been verified perturbatively in the coupling, perturbatively in d — 4 where d 
is the number of dimensions of space, perturbatively in 1/N where N = 4 is the case of 
interest, nonperturbatively in numerical simulations, and by experiment for N = 3. 
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(2) 
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scales, and in particular on arbitrarily long length scales. Long wavelength 
oscillations of the chiral order parameter are the defining feature of physics 
near T c , and our goal later in this talk will be to discuss whether they may 
occur in a heavy ion collison and how they can leave an observable signature. 

If the transition is in the 0(4) universality class, the equilibrium physics 
of the order parameter near T = T c at wavelengths long compared to 1/T is 
classical^ and is described by the Ginzburg-Landau free energy 




Ha\ . (3) 

Here /i 2 and A > are temperature dependent and T c is the temperature at 
which /i 2 = 0. We have introduced an explicit symmetry breaking term Ha 
which tilts the potential, selects a a direction, and gives the Goldstone bosons a 
mass. If the underlying microscopic theory which under renormalization flows 
to (Q) in the infrared were in fact that of a magnet, H would be proportional 
to an externally imposed magnetic field. In QCD, H is proportional to a 
common mass m q for the two light quarks. A nonzero H turns the second order 
phase transition into a smooth crossover. If H is small, the theory (§) yields 
universal, quantitative predictions for the behavior of the order parameter 
((f)) as a function of T for T near T c .|8], |], |10|, |ll|, [12| The critical exponents 



describing physics at H = 0, T ^ T c and at T = T c , H = have been 



computed to high order in a perturbative expansion in A, [13 and by numerical 



simulation. [|14|] The full equation of state — namely (0) as a function of T and 
H — has been obtained numerically by Toussaint. [|T^ 

Finite temperature lattice QCD simulations can test whether the QCD 
phase transition is in fact in the 0(4) universality class. Present simulations 
provide evidence that in two flavor QCD the phase transition is not first order, 
and is plausibly second order in the chiral limit. 0, |j As one example, let us 
consider the results of |J. They compute the order parameter as a function 
of T over a range of temperatures near T c for two values of m q , the lighter 
of which corresponds to m^/nip ~ 0.3. They find that the order parameter 
decreases rapidly but smoothly over a range of temperatures about 10 MeV 
wide centered about a temperature T c ~ 140 — 160 MeV. They then show that 
their results can be fitted equally well using the 0(4) equation of state |l2f and 



the mean field equation of state. However, the two fits make very different 

4 Equilibrium finite temperature quantum field theory is formulated in 4 dimensional 
Euclidean space, where the fourth dimension is periodic with period 1/T. Integrating out 
modes with energy ~ 1/T and higher, yields an effective three dimensional classical field 
theory describing physics at wavelengths longer than 1/T. 
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predictions for behavior at smaller quark masses. Thus, the present simulation 
is consistent with the hypothesis that the transition is in the 0(4) universality 
class, but until smaller quark masses are explored no stringent tests will be 
possible. To see the universal long wavelength physics one must use light 
enough quarks that m n {Tc) < T c , where m 7r (T c ), the pion inverse correlation 
length, is of order or slightly more than m^. In the simulations of 0, m n > T c . 
This means that the correlation length is just not long enough for the physics 
to become effectively three dimensional and classical as in (|3|). In present 
simulations done by a number of groups]^, |(| which have all had m 7T /m p ~ 0.2 
or more, although a variety of methods of extracting critical exponents have 
been explored, it has not yet been possible to cleanly distinguish between, say, 
0(4) and mean field exponents.^] 

We have seen that the long wavelength oscillations of the orientation of the 
condensate are central to the physics of the chiral phase transition. Although 

5 It is worth enumerating the logical possibilities which could lead to a failure of the 
hypothesis that a second order phase transition with 0(4) exponents will be seen in lattice 
simulations. We only know that an infrared fixed point with the appropriate symmetry 
exists. We do not know that QCD is in fact in the basin of attraction of this fixed point. 
Renormalization toward the infrared starting from QCD may lead to the 0(4) fixed point, 
it may miss the fixed point and yield a first order transition, or, least likely, it may lead to 
some as yet undiscovered other fixed point, with different exponents. The first order option 
does in fact arise if the strange quark is lighter than a certain (tricritical||, §,JIo)) value. 
Simulations with two light and one heavier quark done by the Columbia group (lj suggest 
that in nature, the strange quark is heavy enough that no first order transition occurs. It 
is very important to verify this conclusion in the next generation of lattice simulations, 
particularly as it has recently been questioned. |l6) A final logical possibility which would 
lead to nonobservation of 0(4) exponents is that the Ginzburg region may be too small to 
see the true critical behavior. Fluctuations of the order parameter are important only close 
to T c — for \t\ = \T— T c \/T c < £g where tc like T c , is not universal. Outside the Ginzburg 
region, that is for \t\ > tc, mean field theory is valid. If there is an appropriate small 
parameter in the theory, to can be small. For example, Jl7| in ordinary, low temperature, 
BCS superconductors, to ~ (T c /Ef) 4 where the Fermi energy Ef is 10 3 to 10 4 times bigger 
than T c . As another example, Jl^] numerical simulations of the 2+1 dimensional Gross-Neveu 
model suggest that to is small in this theory. Here, the small parameter is 1/N, where N is 
the number of fermions in the theory. In the large N limit in this theory, the phase transition 
remains in the Ising universality class but tc goes to zero like 1/iV.Q (Note that in QCD, 
increasing either N c or Nf has much more drastic effects like making the transition first 
order or changing the symmetry of the order parameter or destroying asymptotic freedom.) 
There is no parameter in QCD, small or not so small, which when taken to zero reduces 
tc, leaving the transition otherwise unaffected. We can therefore expect tc to be of order 
one. To exclude the possibility that it is unexpectedly small, one needs lattice simulations 
with light enough quarks that ~ T c . This should suffice to distinguish between 0(4) 
and mean field exponents. 
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measuring critical exponents, and thus gaining a quantitative understanding 
of this physics, is possible on the lattice, it is extremely unlikely that such 
measurements can be done in heavy ion collision experiments. However, as I 
explain in the rest of this talk, classical long wavelength pion oscillations may 
arise in a heavy ion collision, and can leave a signature.^, |K], |21[ 



A long wavelength oscillation of the order parameter consists of large re- 
gions in which the chiral condensate points in directions other than (a, 0). Let 



us first consider an idealized situation [^, |23|, g4| in which there is a single 
large region in which is uniform in space and misaligned. Because the pion 
mass is nonzero due to the explicit chiral symmetry breaking introduced by 
nonzero quark masses, in such a region of disoriented chiral condensate the 
field would oscillate about the a direction. That is, say, 

= w(cos 9, 0, 0, sin 6) with 6(t) ~ sin(m 7r t) . (4) 

In this idealized case, is independent of x. This can be thought of as an 
infinitely long wavelength oscillation of the orientation of the condensate, un- 
polluted by short wavelength clutter. (In a heavy ion collision, we will see 
that if long wavelength oscillations arise, they occur superposed with short 
wavelength "noise" and so a smooth region of disoriented chiral condensate is 
an idealization. Bjorken et al have suggested that circumstances closer to the 



idealized case may actually arise in hadron hadron collisions. [24]] ) 

An idealized disoriented region in which the disorientation is in the iri — 712 
plane corresponds to all charged pions (equally positive and negative since the 
fields are real), while if the disorientation is in the ir^ direction as in (^) it 
corresponds to all neutral pions with no charged pions. More generally if we 
define the ratio 

R = ^ (5) 

n n o + n n + n - 

then each disoriented region will yield pions with some fixed R. Under the as- 
sumption that all directions of disorientation are equally probable, an ensemble 
of "events" will yield a probability distribution for R [^, [24], |2T 



P(R) = -^. (6) 
2VR 

Note, for example, that the probability that the neutral pion fraction R is less 
than 0.01 is 0.1! In "events" with, say, 100 pions, the distribution (||) is very 
different than the binomial distribution obtained if each individual pion was 
independently randomly neutral or positive or negative. Rather than rolling 
the dice once per pion, and getting an R distribution looking like a Gaussian 
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about R = 1/3, in a population of idealized "events" just described, the dice 
are rolled once per event yielding the much broader, and skewed, distribution 

Now with a hint of how long wavelength oscillations of the order parameter 
could leave a signature, we return to the question of whether they arise. An 
equilibrium second order transition seems ideal, as near the critical tempera- 
ture the correlation length is infinite and oscillations occur on arbitrarily long 
wavelengths. However, because the up and down quarks (and hence the pion) 
are not massless, we have seen that the putative second order transition is 
smoothed out, and the correlation length is not long compared to 1/T C , and 
classical oscillations do not arise. Although we are used to thinking of nature 
as being close to the chiral limit, because m^jm p is small, it is unfortunately 
not close enough for m n /T c to be small. On the lattice, one can envision using 
quark masses smaller than in nature, and so getting closer to having a second 
order transition. In an experiment, one does not have this option. 



Fortunately, this is not the end of the story. [ 3D | There is no reason to as- 
sume that during a heavy ion collision the long wavelength modes of the order 
parameter stay in thermal equilibrium as the plasma expands and cools and 
chiral symmetry breaking occurs, even if local thermal equilibrium at a high 
temperature is achieved at early times. There is a simple non-equilibrium 
mechanism, which could plausibly operate in a heavy ion collision, which nat- 
urally leads to great amplification of long wavelength modes, relative to what 
would be seen in thermal equilibrium. 

Let us consider an idealization that is in some ways opposite to that of 
thermal equilibrium, namely the occurrence of a sudden quench from high to 
low temperatures in which the (a, tt) fields are suddenly removed from contact 
with a high temperature heat bath and subsequently evolve according to zero 
temperature equations of motion.f] Unlike in the equilibrium case where uni- 
versality was our guide, away from equilibrium we must make a non-universal 
choice of model Lagrangian to obtain equations of motion for the order pa- 
rameter. The linear sigma model 

C = J d A x{^rd^ a - - {f(j> a - v 2 ) 2 + Ha\ , (7) 



6 Quenching is an ad hoc assumption, not expected to be realized in a heavy ion collision 
in any detail, and so it is reasonable to consider different dynamical assumptions as has 
been done in (26[ ^tJ. The hope is that qualitative behavior observed after a quench is 
representative of phenomena occurring in realistic conditions in heavy ion collisions in which 
the long wavelength modes are not in equilibrium. 
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is a reasonable choice p0| although other models have also been studied. [ETf In 



20(1 , the nonlinear classical equations of motion derived from this Lagrangian 
were solved numerically on a lattice with spacing a = (200MeV) _1 , and with 
parameters in the potential chosen such that f m and m, have their zero tem- 
perature values and m a = 600MeV. In the initial conditions, (<p) ~ and the 
orientation of <j) was chosen randomly at each lattice site, as appropriate for 
initial conditions above T c at a high enough temperature that the correlation 
length is shorter than a. We chose ((f) 2 ) 1 / 2 = v/2, a choice which has since 
received some justification from the work of Randrup, as we discuss below. 
Because of the explicit symmetry breaking, is soon oscillating about the 
a direction everywhere in space. Upon Fourier transforming, one finds that 
modes of each component of if with wave vector k oscillate about tt = with 



io ~ yk 2 + m 2 , as expected. The behavior of the amplitude of the low mo- 
mentum pion modes as a function of time is striking. The initial conditions 
had a white noise power spectrum, with all modes having equal amplitudes. 
At very late times, things become boring once again, as the system approaches 
an equilibrium configuration in which equipartition of energy holds. Over a 
wide range of intermediate times of order 5 to 50 times m~ 1 , however, the long 



wavelength modes are greatly amplified. (This was confirmed in p8| |29|.) The 
affected modes are those with \k\ less than about m^. The longer the wave- 
length, the bigger the amplification; the longest wavelength mode not affected 
by the finite size of the box has its amplitude squared amplified by more than 
a factor of 1000 relative to that of the short wavelength modes which were 
not amplified at all, and by about a factor of 50 relative to its value at late 
times in equilibrium. For a visual image, think of long wavelength ocean swells 
superposed with lots of short wavelength chop. The essential qualitative phe- 
nomenon whose cause and consequences we now discuss is that long wavelength 
oscillations of the orientation of the condensate (pion oscillations) have been 
excited — the short wavelength modes have not been damped out. One does 
not see smooth domains, and the power spectrum is not well characterized 
by a single correlation length or domain size.|T(| This picture is more compli- 



cated than the idealized "smooth" region of disoriented chiral condensate we 
discussed above. 

The emergence of long wavelength pion oscillations is a striking qualitative 
phenomenon, for which there is a simple qualitative explanation. |J20| One can 
linearize the equations of motion for the pion field to obtain 

d 2 - 

— n(k,t) = -m 2 cfi (k,t)n(k,t) , (8) 
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where the time dependent "mass" is given by 



-Xv 2 + k 2 + X(<p 2 )(t). (9) 



Linearizing in this way constitutes an uncontrolled truncation, but it does 
yield some insight into the behavior found by numerical solution of the non- 
linear equations. One sees that the non-equilibrium dynamics leads to a time 
dependent m 2 a which can be negative for modes with small enough k. This 
means that these modes will experience periods of time during which they are 
unstable to exponential growth. The longer the wavelength, the greater the 
amplification. In an equilibrium phase transition, explicit symmetry breaking 
keeps the correlation length at T c too short to be of interest. This, then, is 
a simple non-equilibrium mechanism which leads to the amplification of arbi- 
trarily long wavelength modes of the pion field even though the pion mass is 
nonzero. 

I will return to analyzing the results of the simulation just described in 
terms of the ratio R, but let me first sketch some improvements. The most 
important effect which has been left out above is expansion. The late time be- 
havior with no expansion (equilibration of the system) is qualitatively different 
than that with expansion. Expansion causes the energy density to drop, and 
the dynamics linearizes as the modes stop interacting at late times. Hence, 
at late times the ratios of amplitudes between different modes stops changing. 
This is the classical analogue of "freeze out". The hope in |20| was that once 



expansion is included, the modes would freeze out with amplitude ratios as 
found at times of order 5 — 50 m' 1 , namely with long wavelength modes am- 
plified. This was in fact found in some,|30|, |31], |32| but not all,|29j subsequent 



simulations which began with initial conditions akin to those above and which 
included effects of expansion in several different ways. (For treatments in other 
contexts which include effects of expansion, see f33j.) 

Another goal in implementing expansion, however, is to relax the quench 
assumption. This has been implemented by Randrup.|34j] He begins with a 
linear sigma model configuration in thermal equilibrium at 400 MeV,[] and 
adds a damping term to the equations of motion which damps energy from 
the system as if the system is undergoing a three dimensional boost invariant 
expansion. Initially, (0) is small, but because of the large initial temperature 



is much larger than in the initial conditions of ||20|| . Randrup finds an 



7 Although this is certainly less ad hoc than the initial conditions in [Q, it is not yet 
realistic. At high temperatures like these, at wavelengths of order 1/T there are presumably 
degrees of freedom other than those of the linear sigma model which are important. 
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initial rapid decrease in (0 2 ), during which ((f)) hardly changes and (0 2 ) de- 
creases to ~ v 2 /2\ Thus, the expansion takes the system to a configuration 
similar to that configuration chosen in 0] as an initial condition for evolution 
after a quench, which in turn was supposed to model the effects of expansion. 
The subsequent evolution in Randrup's simulation is similar to that in PU 



long wavelength modes grow, although because of the expansion they later 
freeze out as described above. Thus, Randrup's result provides some support 
for some of the ad hoc choices made in |2(J . 

A more ambitious project than the classical simulations described so far is 
to include the effects of quantum fluctuations (in addition to classical thermal 
fluctuations) on the dynamics of the order parameter. We have seen that with 
the pion mass nature has dealt us, the equilibrium correlation length is not 
long compared to 1/T C , and this suggests that quantum effects may matter. 
However, the nice thing is that the instability discussed above is an instability 
towards growth of long wavelength modes, which as they grow in amplitude 
become more classical. In equilibrium, the occupation number of a mode 
is ~ T/uj which is ~ T/m^ at long wavelengths, so any growth relative to 
equilibrium pushes the mode into the classical regime. Simulations including 
quantum fluctuations have now been done by several groups, |]3T], and 
they in fact find the same instability to growth of long wavelength modes 
found classically. However, in quantum treatments done to this point it has 
been necessary to make linearizing approximations. I argue in |T(| that the 
simulations of [[21J show that the growth of long wavelength modes in the full 
nonlinear theory is in fact greater than in the linearized theory. Hence, it is best 
to view the two treatments as complementary. The quantum treatment can 
be used to study under what circumstances the instability is present, but once 
the instability sets in and long wavelength modes begin to grow, a classical 
treatment becomes appropriate. Happily, the two approaches seem to yield 
approximately similar outcomes. 

Because a non-equilibrium transition is required, the theoretical founda- 
tions of all the simulations are shaky. Using the linear sigma model to de- 
scribe the dynamics seems reasonable for long wavelength modes, but cannot 
be justified in a controlled way. Either neglecting quantum effects or includ- 
ing them in a Hartree approximation risks missing some of the physics. Most 
important of all, we do not know what initial conditions to impose on the 
long wavelength modes of the chiral order parameter at the time when par- 
tonic language ceases to be appropriate and a description in terms of the linear 
sigma model becomes appropriate. Nobody can do a simulation which follows 
non-equilibrium dynamics beginning with cascading partons and ending with 
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long wavelength oscillations of the condensate. This makes it impossible to 
predict that long wavelength pion oscillations will be amplified in a particular 
experimental setting. The way to know that this phenomenon has occurred is 
to detect it, and to this we now turn. 

The defining signature^ of the presence of long wavelength classical pion 
oscillations is large fluctuations in the ratio R. If a particular k mode ends up 
oscillating in the 7r 3 direction, then it will become neutral pions moving in the 
direction of k. Similarly, oscillations in the n 1 — n 2 plane will become charged 
pions. If this were the only physics occurring, one could bin the data as a 
function of pseudorapidity 77 and azimuth ip and the bin by bin distribution 
of R would be (||). The problem, of course, is the large background due to all 
the other pions in the event. 

What is needed is a detector which covers as much of 77 and ip as possible, 
which counts both charged pions and photons, which accepts pions down to 
as low Pt as possible (50 MeV or even lower; this may mean running with 
reduced magnetic field), and which accepts photons down to pt ~ or so. 

The detector should be segmented in 77 and <p finely enough that the number 
of photons and charged pions in each segment can be accurately counted. 
Where should one be in 77? Central rapidity is cleanest for a theorist, but 
at RHIC low pt photons are hard to count there, so perhaps somewhat off 
central rapidity may be better. Accurate measurement of px and particle 
identification is not necessary, since at RHIC most charged particles are pions 
and most photons are from 7r 's. However, it would be very helpful to have 
enough p? information to make even a crude cut, keeping only pions with low 
Pt, say pt < 2 — 3 m^. The data, then, consists of an event by event catalogue 
of the positions in (77, ip) of the photons and charged pions. From this, one can 
compute R = n 7 / (n 7 + 2n c h) in bins in (77, <p). 

In thinking about how to analyze real data, it is an instructive exercise |pE0| 



to go back to the simulations of [g(J and make explicit how the amplification 
of long wavelength modes is manifest in the number ratio R. For each of the 
three 7?, and for each k on the momentum space lattice with spacing 2-K/a, 
define the "particle number" 

nUt) = - iun\k, t) + r(k, t) 2 , (10) 



8 Other signatures which have been discussed include an excess in the total number of 
pions at low pr independent of their charge, |3lJ [32[ |36|. |37[| various effects on pion pair 
correlations, |58] and various electromagnetic effects. |39|, |37| Any of these would be telling 
if seen in conjunction with large fluctuations in R, but it is not clear that they would be 
unambiguous if seen alone. 
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Fig. 1. R as a function of angle in momentum space, binned in 10° x 10° 
bins. The top panel is at t = 160 fm when nothing interesting is going on, and 
illustrates the magnitude of "ordinary" fluctuations. The bottom panel is at 
t = 22 fm, a time when long wavelength modes have been amplified. 

where ir l (k,t) and n l (k,t) are the Fourier transforms of Tr l (x,t) and its time 
derivative. If we had included expansion, had waited until the system lin- 
earized, and if we assumed that the classical fields correspond to quantum 
mechanical coherent states, these would be the mean particle numbers per 
mode at late time. Although the n's should really be viewed just as time de- 
pendent quantities defined by (|T0|), I will refer to them as numbers henceforth. 
The lattice has 64 3 points, so we now have the number of neutral pions and the 
total number of pions at 64 3 points in momentum space as a function of time. 
I now make a cut, keeping only those modes with \k\ < 300 MeV. For each 
remaining k, I compute the polar angle 9 and the azimuthal angle <p describing 
the direction in momentum space. (In an experiment, one would use rj. Here, 
without expansion, 9 is the right variable.) I then divide momentum space 
into 10° x 10° bins in (0, ip), and discard the bins within 30° of the north and 
south poles, as they each contain too few modes. For each of the remaining 

— * 

bins, I add up the number of neutral pions in all the modes with fc's in the 
bin, do the same for the total number of pions, and compute R. 

The results are shown in Figure 1, at two different times. The top panel is 
at a very late time, when the system has equilibrated and nothing interesting 
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Fig. 2. Histograms of the distributions of R in Figure 1. The shaded histogram 
corresponds to the lower panel in Figure 1, and the histogram drawn simply 
ELS £1 line corresponds to the upper panel. 

is going on. The bottom panel is at an intermediate time when the long 
wavelength modes have large amplitudes. It is quite clear in the figures that 
the R distribution is much broader in the lower panel. There are more high 
spots and low spots. Figure 2 shows histograms of the distribution of R values 
in the 432 bins in Figure 1. The unshaded histogram corresponds to the top 
panel, and the shaded histogram corresponds to the bottom panel, in which 
long wavelength oscillations are present. We see that the effect of the long 
wavelength modes is to broaden the distribution (look at the tails), and skew 
it to the left (look at the peaks). This is exactly what one would expect from 
an admixture of a 1/yfR distribution. Note that if the 300 MeV cut on \k\ 
is pushed too high, the short wavelength modes swamp the long wavelength 
modes and both distributions look the same. If the cut is too low, there are 
very few modes per bin, and both distributions become wider. Eventually, if 
the cut is too low (or if the bin size is too small) and there is only one mode per 
bin, then both distributions are l/\/R. Increasing the bin size above 10° x 10° 
narrows both distributions by the same factor. The central limit theorem 
"acts" on both distributions as the bin size is increased, but the presence of 
large amplitude long wavelength modes "delays its action" . |41] The choice of 
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cuts and bin sizes in Figs. 1 and 2 was driven by the background (short 
wavelength modes) and statistical fluctuations (arising if there are too few 



modes per bin) in the simulation in [20]. In analyzing real data in an analogous 



fashion, such choices will be driven by the background (high pions) and 
statistical fluctuations (arising if there are too few pions and photons per bin) 
and the fact that one counts photons rather than 7r°'s. Hence, the choices I 
have found convenient need not be the ones which are convenient in analyzing 
real data. 

With real data, one should try bins of all sizes in both 77 and ip, ranging from 
bins which on average contain a few tens of pions, to bins in rapidity covering 
27r in ip, to "bins" which consist of whole events, and construct R distributions 
as above for each bin size. (The data should also be analyzed using the wavelet 
formalism of pT|.) It is also very important to look for singular bins far out 



on the tails of the various distributions. The simulations of [|3U| suggest that 
events of this sort may occur, although they are presumably rare. Broadened 
(and perhaps skewed) R distributions may be more common. Both should 
be looked for as functions of increasing beam energy and projectile size, and 
decreasing impact parameter. 

To what should one compare an R distribution obtained from data? One 
can use an R distribution obtained from an event generator in which pions 
are independently randomly neutral, positive, or negative. If there is no sig- 
nal in the data, this should fit very well, because any badness of fit in each 
individual multiplicity distribution should cancel in the ratio R. However, it 
may be preferable to do the analysis without reference to an event generator. 
Steinberg has suggested creating a data set by rotating the charged pions 
in each event in azimuth relative to the photons, say by 90°. One would look 
for broadening and skewing of the histograms from the real events relative to 
those from the rotated events. 

Although there are fewer complications at RHIC than at lower energies, 
given the lack of theoretical certainty it is certainly wise to look today. WA98 
at CERN is in the process of analyzing their event by event charged particle and 
photon data. [53] They are looking for, and so far have not found, PSI singular 
events which behave as if more than ~ 30% of all the pions in the event come 
from a single idealized region of disoriented chiral condensate. They are in the 
process of binning the data in (r/, ip) and looking for fluctuations in R using 
various bin sizes. Even before they finish their analyses, we can learn useful 
lessons for future experiments. Being able to cut on pt, and get rid of the 
high momentum pions, would really help, as would having a charged particle 
veto in front of their photon multiplicity detector. NA49 is planning to search 



13 



for fluctuations in the number of charged pions only, as they do not count 
photons. mf] PHENIX, PHOBOS, and STAR all plan to search for unusual 
fluctuations in R at RH1C. 

The conclusions from theory for experiment are qualitative. Consider a 
heavy ion collision which is energetic enough that there is a central rapidity 
region of high energy density and low baryon number. If, as this region cools 
through the chiral transition, long wavelength modes do not stay in equilib- 
rium, then there is a robust mechanism which can lead to the amplification of 
long wavelength pion oscillations. This cannot happen in thermal equilibrium. 
One should not take the details of any of the theoretical simulations too seri- 
ously. All make dynamical assumptions and all rely on guesses for their initial 
conditions, since nobody can go from initial parton dynamics to the dynamics 
of the chiral order parameter at long wavelengths. What one should take seri- 
ously is that the place to look for signatures of the dynamics of the chiral order 
parameter during the phase transition is the low pt pions. We have discussed 
one effect, leading to such a signature, but perhaps the best outcome would 
be the detection of completely unexpected phenomena at low p?. Long wave- 
length oscillations of the condensate would be observed by detecting unusual 
fluctuations in R, the number ratio of neutral pions, event by event as a func- 
tion of rapidity and azimuth. One should look both for broadened and skewed 
R distributions, and for singular events with particularly large fluctuations. 
Experimentalists are trying to disorder the chiral condensate and then detect 
its disoriented oscillations, as it returns toward its ordered ground state. They 
will continue at RHIC and the LHC, where sufficient energy should be avail- 
able. Detection of unusual fluctuations in R in a heavy ion collision would be 
a dramatic and definitive signature of an out of equilibrium chiral transition. 
The ball is in the experimentalists' court, and we wish them well. 
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9 MINIMAX is looking for regions of disorientated chiral condensate at high rapidity 
in hadron-hadron collisions at Fermilab. The dynamics in these events is surely different 
from that we are discussing, but the signatures are similar and therefore their analysis 
techniques should prove useful when generalized for use in larger acceptance experiments. 
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